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Using the hadronic matrix elements from the lattice, Bk and £ s , involving only the 4-quark 
operators for A flavor = 2 Hamiltonian relevant for K — K, Bd — Bd and B 3 — B 3 mixing, along with 
Vcb, we deduce a non-trivial constraint on the SM, sin(2/3) = 0.87 ± 0.09. This deviates from direct 
experimental measurements via the tree process, b — > ccs as well as the one via the penguin-loop 
b — » s decays by around 2.1 and 2.7 a respectively. If these deviations are confirmed they would 
imply the presence of new physics rather pervasively in both Bd — Bd (i.e very likely in b — > d) as 
well as in b — > s transitions requiring a beyond the SM CP-odd phase. Consequently, improvements 
in the relevant lattice calculations should be given a high priority. 



In the past several years spectacular performance of 
the B-factories has led to a milestone in our understand- 
ing of CP violation phenomena. Results from the B- 
factories, along with Standard Model (SM) predictions, 
based in large part on lattice calculations of relevant weak 
matrix elements, provided a striking confirmation of the 
Cabibbo-Kobayashi-Maskawa (CKM) [J [3 paradigm [3|. 
For the first time it was experimentally established that 
the single (CP-odd) phase of the KM ansatz is able to 
account to an accuracy of about 15% the observed 0(1) 
CP asymmetry seen in the gold-plated, B — > iJjKs type 
of modes as well as the minuscule (i.e. O(10 -3 )) CP 
asymmetry seen in Kl decays long ago [J]! While this 
is a remarkable success of the CKM picture, we should 
note that there are essentially compelling theoretical ex- 
pectations that beyond the Standard Model (BSM) CP- 
odd phase(s) must also exist in nature. Finding evidence 
for these is clearly an important challenge for Particle 
Physics. 

One strategy to search for these new phase(s) is to 
look for deviations from the SM in its precise descrip- 
tion of the asymmetry as seen in (say) B — > ipK s decays. 
Since in this case one is clearly looking for effects that 
are small in comparison to the dominant contribution of 
the SM, precision in theory and in experiment becomes 
highly desirable. Another class of searches are amongst 
those channels wherein the SM predicts vanishing asym- 
metries; these belong to a class of null tests [J [B| for 
the SM. Sizable CP asymmetries in such channels would 
then have to be ascribed to BSM sources. In this work, 
we will show that using clean input from lattice calcula- 
tions leads to extremely interesting indications of beyond 
the SM CP-odd phases (s) in both of the aforementioned 
class of tests. Furthermore, there are good reasons to 
expect that these lattice calculations can be improved 
further to provide even more stringent tests. 

In this work, using primarily the input from the lat- 
tice for the non-perturbative matrix elements of the four- 
quark operator of the A flavor (F) = 2, SM Hamiltonian, 
we show that the value of sin (2/3) that we thus determine, 



with this clean input, exhibits some deviations from that 
measured directly in B-factory experiments. These dis- 
crepancies are potentially of crucial importance as they 
signify the possible presence of physics beyond the Stan- 
dard Model in Bd — Bd mixing, as well as in b — > s pen- 
guin transitions. These observations lead to the impor- 
tant conclusion that the effects of physics beyond the SM 
are possibly making their presence felt rather pervasively 
both in b —> d as well as in b — > s short distance penguin 
transitions. These deviations are specially important as 
they are obtained without making any use of V u b, which 
has been very problematic in recent years as the value de- 
duced from inclusive methods disagrees appreciably with 
that deduced from the exclusive approaches Q . 

To reiterate, we show that, along with V cb , only the 
hadronic matrix elements of one type of 4-quark opera- 
tor: [/i7 M (l — 7s)<z] 2 , with h = s (q = d) and h — b (q = s 
and d) are needed for providing a non-trivial constraint 
on the value of sin(2/3) in the SM. Note that for all three 
of the relevant mixing-matrix elements, for K, Bd and B s 
mesons, in fact the strange and the 6-quarks can be con- 
sidered very heavy compared to the light u and d quarks, 
explicitly demonstrated recently by the RBC-UKQCD 
collaborations' determination of Bk |7j. Thereby only 
SU{2) x SU{2) (rather than ST/ (3) x SZ7(3)) chiral pertur- 
bation theory need be used for chiral extrapolations from 
the light quarks on the lattice to their physical value re- 
sulting in significantly improved accuracy. Furthermore, 
it is important to note that from B-mixings in fact we 
only use the SU (3) breaking ratio of mass differences, 
AMbJAM^ [8(. The prognosis for steady improve- 
ment in the accuracy of these matrix elements therefore 
seems to be quite good. 

We now present the analysis based on the thinking 
outlined above. Fig.[T]shows the region of the (p, fj) plane 
allowed by |F c6 |, AM Bg /AM Bd and e K . No use of V ub 
is made here; this is an important difference from our 
previous work Q . This is done to alleviate any concerns 
that the problems afflicting V u b that we alluded 

to above, may be unduly effecting the results. The key 
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FIG. 1: Unitarity triangle fit in the SM. All constraints are 
imposed at the 68% C.L.. The solid contour is obtained us- 
ing the constraints from ek and AA/s s /AA/s d . The regions 
allowed by a^K and a^ +n i + 2K s )K s are superimposed. 



ek = (2.232 ±0.007) x fO" 3 

Am Bs = (17.77 ±0.10 ±0.07) ps" 1 [22] 

Am Bd = (0.507 ±0.005) ps" 1 

\Vcb\ = (40.8 ±0.6) x 10~ 3 

B K = 0.720 ± 0.013 ± 0.037 [7] 

£s = 1.20 ± 0.06 

A = 0.2255 ± 0.0007 [23] 

m t , poie = (170.9 ± 1.8) GeV [24] 

m c (m c ) = (1.224 ±0.017 ±0.054) GeV [25J 

771 = 1.51 ±0.24 [26] 

772 = 0.5765 ± 0.0065 [27] 

r) 3 = 0.47 ± 0.04 [28] 

TABLE I: Inputs that we use in the unitarity triangle fit. 
When not explicitly stated, we take the inputs from the Par- 
ticle Data Group Q. 



inputs being used are as follows: 

• The K — K mixing matrix element conventionally 
parametrized as Bk 0] needed from the lattice to 
make use of the indirect CP-violation parameter, 
ex from Kl — > 7T7t, 



B 



K 



0.720 ±0.013 ±0.037 



(1) 



• The SU(3) breaking ratio, £ s , from the lattice, 
needed in conjunction with Ams s / ArriB d , which 
is the ratio of the B s and Bd mass differences 11 
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h 



'Bj, 



= 1.20 ±0.06 



f Bd VB d 



(2) 



• As far as V c b, which is another input that is neces- 
sary, we note that there is a ~ 2a tension between 
the extraction of \V c b\ from inclusive and exclusive 
decays 14. 15j|: 



41.7 ±0.4 ±0.6 inclusive 



r 41.7 

1 38.7 



ik,ix io* = < — — - — ; ( 3) 



Therefore, in the numerics we use the weighted av- 
erage of these two determinations: 



\V cb \ = (40.8 ±0.6) x 10" 



(4) 



Thus from Fig. Q] it is clear that even without the in- 
clusion of \V u b\, the prediction for sin(2/3) deviates from 
the experimental determinations summarized in Table HT1 
From the chi-squared minimization we find: 



sin(2/3)]^ c b tion = 0.87 ± 0.09 



(5) 



Our fitting procedure consists in writing a chi-squared 
that includes all experimental measurements and lattice 



determinations; this implies that we assume gaussian er- 
rors. The input values that we use in the fits are summa- 
rized in Table |H The SM expressions for AM B J AM Bd 
and ek can be found, for instance, in Ref. [16[. Follow- 
ing Ref. [l3|, we include in ek the term proportional 
to the / — component of the K — > tttt amplitude, 
whose contribution is effectively taken into account by 
the multiplicative factor K e . The calculation of the lat- 
ter is affected by non-perturbative uncertainties and, fol- 
lowing the analysis in Ref. 17, [HI, [l|| [2(|, we take 
k £ = 0.92 ±0.02 [Hj]. 

The constraint in Eq. |5J deviates by 2.1c from sin(2/3) 
extracted from the tree-level decay B — > J/ipKs and by 
2.7a from the average of the three penguin-dominated 
B -► ((j),ri',KsKs)Ks modes (see Table HU). We have 
chosen to concentrate on the time-dependent CP in 
these three b — > s penguin-dominated modes (amongst 
many) as QCD factorization as well as several other ap- 
proaches show that for these the QCD corrections are 
very small [2£" 
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|32j|. In passing, we do want to 
draw attention, though, to the fact that there are many 
more b — > s penguin-dominated modes 33], and curi- 
ously in most of these modes, the central value of the 
time dependent asymmetry is below the predicted SM 
value of Eq. (O [9( or for that matter the value directly 
measured via B — ► tpK s , 0.681 ± 0.025, indicating again 
possible problem with the SM description of the b — * s 
penguin dominated modes. Thus the fact that these data 
show indications for the need for a BSM-CP-odd phase 
reported earlier Q does not appear to originate from the 
use of a faulty input for V u b- 

Note that in part this clean constraint Eq. ([5]) has 
been made possible without any use of V u b due to the 
significant (almost a factor of three) reduction in errors 
in lattice determination of Bk achieved in the past few 
years 34] . In Table |TT] we compare the prediction |5]| to 
the various tree and penguin modes. While these penguin 
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FIG. 2: Unitarity triangle fit in the SM. All constraints are 
imposed at the 68% C.L.. The solid contour is obtained using 
the constraints from ek, AMb 3 / AMs d and |Vub/Vcb|. The 
dashed contour shows the effect of excluding |Kib/V c b| from 
the fit. The regions allowed by a^K and a^ r j, +v i +2 K s )K s are 
superimposed. 



modes have been already getting considerable attention 
in the past « 2 years, our analysis makes clear that the 
low value of sin(2/3) that they yield compared to the SM 
constraint in Eq. has nothing to do with the difficul- 
ties affecting V U b- Furthermore, we find that the difficul- 
ties of the SM description of CP violating B-decays may 
not just be confined to the b —> s penguin modes but in 
fact may even be there in the "gold-plated" B — > tpK s 
mode. Since the latter is a combination of Bd — Bd mix- 
ing and the tree decay, b — > ccs, it is natural to be more 
suspicious that BSM physics in A6 = 2, loop process, 
may be causing this deviation in the former, though the 
tree decay may also be in part the cause. 

We now turn our attention to V U b to see how it can 
change the above results. As has recently been empha- 
sized by Neubert [io| . the extraction of \V u b\ from inclu- 
sive decays is extremely sensitive to the precise determi- 
nation of the b quark mass. Here we follow Neubert 's 
analysis and do not use B — ► X s j for the extraction of 
TOft . The most recent determinations of | V u b | , from inclu- 
sive and exclusive works are [l(J 3^, 3(|: 



37.0 ±1.5 ±2.8 inclusive 
\V ub \ x 10 4 ={ 35.5 ±2.5 ±5.0 HPQCD (6) 

37.8 ±3.0 ±3.4 ±2.5 FNAL/MfLC 

The two determinations of \V u b\ from exclusive decays 
are based on the same experimental data and we chose 
to adopt the Fermilab/MILC [3F ' 
the semileptonic B — > 7r form factor. 



37j determination of 
Combining the 



inclusive and exclusive results we obtain: 
\V ub \ = (37.2 ±2.7) x 10" 4 . 



(7) 



which is now in good agreement (within 0.7er) with all 
three central values given above. 

Combining that with the inclusive and exclusive ex- 
tractions of V cb we obtain: 



mode 


experiment 


no Vub 


with Vub 




0.681 ±0.025 


2.1 a 


1.7 a 


a.4,K s 


0.39 ±0.17 


2.5 cr 


2.1 a 


a v 'K s 


0.61 ±0.07 


2.3 cr 


1.8 cr 


aK s K s K s 


0.58 ±0.20 


1.4 cr 


0.9 a 


a(4, +v '+K s K s )K s 


0.58 ±0.06 


2.7 cr 


2.5 a 


a(i,+4>+ri'+K s K s )K s 


0.66 ± 0.024 


2.3 cr 


2.1 a 



TABLE II: Experimental determinations of sin(2/3) in b — > 
ccs and b — > sss decays [3^]. Also shown are the deviations 
from the SM prediction obtained without (Eq. (0), and with, 
(Eq. ([9])) the inclusion of V u b in the fit. 



We now present the result we obtain in the full fit (i.e. 
including |V^&|) in Fig.[2j The prediction we obtain reads: 



[sin(2/?)] 



prediction 
full fit 



0.75 ±0.04 



which is consistent, within errors, with Eq. 
with the analyses in Refs. 38Ll33|. 



(9) 

as well as 



Fig. [2] clearly emphasizes a very important role of V ub 
in that it is very useful in excluding a second solution for 
the fit to sm(2(3) in addition to potentially reducing the 
errors appreciably (see also Ref. [!(| for a discussion of 
the second solution). 

The deviation from the experimental results are sum- 
marized in Table ITT1 Notice that no significant change in 
the deviations in B — ► ipK s and in B — > (</>, r/, K S K S )K S 
has taken place from the fit that did not use V u b- 

It is perhaps of some use to extract the values of Bk , 
£ s and V c b that are required to reduce to the 1-cr level the 
discrepancy between the prediction given in Eq. (0 and 
a ('4>+<t>+ri'+KsK s )K s = 0.66 ±0.024. We find that one has 
to choose either B™ w = 0.96 ± 0.04, ^ lcw = 1.37 ± 0.06 
or V cb = (44.3 ±0.6) x 10~ 3 . 

To summarize, in the SM picture of flavor and CP vi- 
olation there seem to be some inconsistencies that might 
be an important hint of new physics at the electroweak 
scale. In particular, the predicted value of sin(2/3) in the 
SM seems to indicate possible inconsistencies with the 
value directly measured via the gold-plated B — > tpK s 
mode and also by the penguin-dominated modes, such 
as B — > (4>,rj' , KsKs)Ks- Furthermore, recall also that 
it seems rather difficult to reconcile the observed differ- 
SD (AAgfS = (14.4 ±2.9) x 10~ 3 ) in the direct CP 



ence 



asymmetries of B° 
prediction AAcp = 



4 K--K+ and B~ 
(2.5 ± 1.4) x 10" 



> K 7r with our 
[H based rather 



\V ub /V c 



cb 



0.0924 ±0.0071 . 



(8) 



loosely on QCD factorization [42J and obtained by allow- 
ing several input parameters to simultaneously take the 
values at the edge of their range. If these difficulties per- 
sist then they would require the need for a beyond the 
SM CP-odd phase. It does not seem like the problems in 
the determination of V u b are the cause of these discrep- 
ancies; so our earlier conclusions regarding the need for a 
new-phase are substantiated 0] ■ We also note in passing 
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the interesting recent work 43J which indicates presence 
of new physics in their analysis of B s — B s mixing (See 
also Ref. [3]). From the perspective of our studies in 
and the current work, which suggest the need for a new 
phase in b — * s, the findings of [43| appear naturally 
related. Therefore, there is a heightened need for fur- 
ther clarifications on these important issues. Improved 
measurements of sin(2/3) via B — > ipK s or via the pen- 
guin modes, e.g. B — > r)'K s , may well have to await 
new experimental facilities, such as a Super-B (or Super- 
Flavor) Factory 45, 4f|. The predicted value of sin(2/3) 
relies heavily on lattice determinations of V c b, the SU(3) 
breaking ratio £ s , the kaon parameter Bk, and V u b- The 
importance of further improvements in these calculations 
can hardly be over emphasized. 
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